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of  the  time-frequency  plane  is  examined.  The  degree  to  which  a  particular  signal  is 
concentrated  is  measured  by  integrating  the  Wigner  distribution  over  the  given  region. 
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associated  operator  is  studied.  Estimates  for  the  eigenvalue  decay  and  the  smoothness 
and  decay  of  the  eigenfunctions  are  established. 
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SECTION  1 
INTRODUCTION 


It  is  well  known  that  the  time-frequency  characteristics  of  a  square  integrable  signal 
cannot  be  arbitrary.  For  example,  no  such  signal  can  be  both  time  and  band  limited. 
The  Heisenberg  uncertainty  principle  provides  another  quantitative  restriction  on  the 
joint  time- frequency  behavior  of  a  square  integrable  signal.  These  facts  indicate  that  a 
signal  cannot  have  all  its  energy  concentrated  in  a  finite  region  of  the  time-frequency 
plane. 

Nevertheless,  in  many  applications  it  is  important  to  use  signals  whose  time- 
frequency  characteristics  are  highly  localized.  Among  other  things,  the  work  of  Landau, 
Pollack,  and  Slepian  [1,2]  produced  a  rigorous  development  of  band  limited  signals  that 
are  as  concentrated  as  possible  within  a  prescribed  timespan.  More  recently,  there  has 
been  interest  in  finding  signals  that  are  localized  in  general  regions  of  phase  space  via 
methods  that  keep  time  and  frequency  on  an  equal  footing  (see  the  papers  of  Daubechies 
and  Paul  [3,4]). 

In  this  paper,  we  study  a  localization  technique  that  uses  the  Wigner  distribution 
to  measure  the  degree  to  which  a  signal  is  concentrated  in  a  particular  region.  This 
leads  to  a  self-adjoint  localization  operator  that  is  easy  to  study  in  terms  of  the  Weyl 
correspondence.  Under  the  Weyl  correspondence,  a  function  of  two  variables — called 
the  symbol — is  associated  with  an  operator  on  functions  of  one  variable.  The  symbol 
of  the  localization  operator  is  simply  the  indicator  function  of  the  given  region  in  the 
time-frequency  space.  The  eigenfunctions  of  this  operator  with  large  eigenvalues  span 
a  subspace  that  can  be  used  to  determine  the  component  of  a  general  signal  that 
is  concentrated  within  the  given  region  of  the  time-frequency  plane:  one  computes 
the  projection  (either  orthogonal  or  weighted  by  the  eigenvalues)  of  a  general  signal 
into  this  subspace.  This  procedure  was  put  forward  by  Flandrin  [5],  who  derived  a 
number  of  useful  results,  including  Lemma  4  below.  It  has  since  been  developed  in 
the  context  of  time-varying  filtering  (see  for  example  the  papers  of  Hlawatsch,  Kozek, 
and  Krattenthaler  [6,7]).  This  paper  is  devoted  to  the  further  study  of  the  localization 
operator  described  above.  The  asymptotic  properties  of  the  eigenvalues  are  studied  and 
it  is  established  that  they  are  C7(A:~3/4).  In  addition,  the  eigenfunctions  with  nonzero 
eigenvalues  are  shown  to  have  faster  than  exponential  decay  in  both  the  time  and 
frequency  domains.  This,  of  course,  leads  to  a  statement  regarding  the  smoothness  of 
these  eigenfunctions.  In  particular,  we  show  that  they  are  analytic.  In  the  last  section, 
some  numerical  examples  are  provided. 
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SECTION  2 

THE  WEYL  CORRESPONDENCE 

The  basic  properties  of  the  Weyl  correspondence  that  we  will  need  are  collected  in 
this  section;  for  a  thorough  treatment,  the  reader  is  referred  to  the  book  of  Folland  [8]. 

Let  /, g  £  L2( R).  A  general  time-frequency  shift  of  /  is 

The  cross- ambiguity  function  of  /  and  g  is 

Mf,9)(T,<r)  =  (p(r,<r)f,g) 

=  J  e*'T<J  e2*iai  f{t  +  r)g(t)dt 

=  J  e2xi<Tsf(s  +  r/2)g(s  -  r/2 )ds. 

The  value  of  A(f,g)  at  a  particular  point  is  the  cross-correlation  between  a  particu¬ 
lar  time-frequency  shift  of  /  against  g.  The  ambiguity  function  is  therefore  a  time- 
frequency  cross-correlation  between  two  functions.  The  Wigner  distribution  is  the 
two-dimensional  Fourier  transform  of  the  cross  ambiguity  function,  thus  giving  it  the 
interpretation  of  a  time-varying  spectrum.  The  Wigner  distribution  can  be  written  as 

W( /,«)((,()  =  J  t-**T(nt  +  T/2)g(t-r/2)dT. 

Several  useful  properties  of  the  Wigner  distribution  are  catalogued  below. 

Theorem  1.  Let  f,g  e  L2( r).  Then 

(i)  W(f,g)(i,  t )  €  i2(R2)  and  \\WU,g)f  =  ll/l|2M2, 

(ii)  W(f,9 )  €  Co(R2)  an d  ||»'(/,s)||0o  <  ||/||||9||, 

(iii)  W(g,f)  =  W(f,g), 

(iv)  W(f,g)(i,t)  =  W(/,j)(«,-().  and 

(v)  W(p(a,b)f,g(a,b)g)U,t)  =  W(f,g)((  -  b,t  +  a). 

The  Weyl  correspondence  uses  the  Wigner  distribution  to  define  a  correspondence 
between  functions  of  two  variables  and  operators  on  L2( R).  It  is  defined,  via  duality, 

by 

(Lsf,g)=  jj  S((,t)WU,gM,t)d(dt 

where  /, g  6  L2{ R)  and  is  a  function  with  appropriate  decay  properties.  S({,t) 

is  the  symbol  of  the  operator  Ls ■  The  following  theorem  of  Pool  [9]  is  useful. 
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THEOREM  2.  A  symbol  S(£,t)  G  L2( R2)  gives  rise  to  an  operator  Ls  that  is  Hilbert- 
Schmidt  on  L2( R).  Moreover,  the  mapping  S  i— ►  Ls  is  a  unitary  operator  from  L2( R2) 
to  the  Hilbert-Schmidt  operators  on  L2( R). 

LEMMA  3.  Ls  is  self-adjoint  if  S(£,t)  is  real  valued. 

The  following  result  was  derived  by  P.  Flandrin  in  [5],  based  on  results  of  Janssen 
[10].  We  will  follow  the  development  in  [8]. 


LEMMA  4.  The  eigenfunctions  of  the  operator  Ls  corresponding  to  a  radially  symmetric 
symbol  S  are: 


21/4  —  1  2  A? 


— 2jti 


PROOF:  Theorem  1.105  in  [7]  shows  that 


W(hj,  hk){(,t)  = 


2(-\)k 


2(  — 1)J 


(2yfizy~k  L[]~k)(U\z\2), 

-2*W\2^z)k->L[k-]){U\z\2), 


for  j  >  k 


for  k  >  j 


where  z  =  t  +  and  L is  the  associated  Laguerre  polynomial.  Set  r  =  |c|  and  note 
that 


( Lshj,hk)  -  JJ  S(r)\V(hj.hk)d(dt 


0. 


( —  1  )J4 7r  /  S(r)e~2rr2 L^iUr^rdr, 


for  j  f  k 
for  j  =  k. 


Hence 

rOO 

A  j  =  ( —  1  )7  4  rr  /  S(r)e~2rr  L*°*(47rr2)  rdr. 

Jo  3 


(1) 
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SECTION  3 

LOCALIZATION  VIA  A  CUT-OFF 


The  localization  operator  we  are  concerned  with  is  Lxn ,  where  xn  is  the  char¬ 
acteristic  function  of  some  bounded  domain  in  the  t  —  £  plane.  We  will  assume  that 
ft  C  [— B,  B]  x  [— T,  T\.  Note  that 

(LxJ,g)  =  J J^W(f,g). 

LXn  is  self  adjoint  since  xn  is  real  valued.  Pool’s  theorem  implies  that  LXn  is  Hilbert- 
Schmidt.  Hence  there  is  an  orthonormal  basis  <f> of  L2(R)  and  real  numbers 
Ai,A2,...  such  that  LXn4>k  =  Xk^k-  The  Hilbert-Schmidt  norm  of  LXn  is  ]T](Afc|2  = 
|n|,  the  measure  of  ft.  We  will  assume  that  the  eigenvalues  are  arranged  in  order 
of  decreasing  absolute  value.  It  is  easy  to  check  that  the  largest  positive  eigenvalue 
corresponds  to  the  maximum  energy  an  L2  function  can  have  within  the  domain  ft.  The 
corresponding  eigenfunction  would  then  be  a  time  function  with  energy  as  concentrated 
as  possible  within  ft.  Our  principal  aim  in  this  paper  is  to  study  the  decay  of  the 
eigenvalues  of  LXn  and  the  smoothness  properties  of  the  eigenfunctions. 

Several  properties  of  the  associated  kernel  will  be  of  importance. 

LEMMA  5.  The  kernel  of  the  operator  LXn  is  given  by  the  equation 

K{s,t)  =  Jx  n(«,i1p)e2l,(,-‘l£^ 

and  has  the  properties 

(i)  K(s,t)  =  0  if  |s  +  t|  >  2T  and 

(ii)  if  the  cross-sections  ft<  of  ft  in  the  £  direction  consist  of  at  most  M 
intervals,  then  |/t'(s,t)j  < 

PROOF:  The  formula  for  the  kernel  is  well  known  (see  [8]).  The  kernel  can  be  written 
as 

where 

F(VJ)=  ( 

Jn, 
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Item  1  follows  from  the  observation  that  Q*  =  0  if  |t|  >  2 T.  We  now  verify  item  2. 
By  assumption,  the  cross-section  is  the  union  of  at  most  M  disjoint  intervals:  fl*  = 
U^Lj £*i>  A]  for  some  integer  L  <  M .  t  herefore  one  can  estimate  that 


1^,01  = 


< 

< 


E 

E 

2  L 


e*Hak+0k)n  sm(n(0k  -  Otk'n) 
nri 

sin(7r (/?*  -  ak)rj) 


< 


Tcrj 

2M 


+  i  tol  +  i 


since  sin(Ar)/x  <  my~j+ •  The  estimate  only  needs  to  be  verified  when  |s  +  t|  <  2T. 
In  this  case, 

\s  -  t\  +  2T  >  \s  -  t\  +  \s  +  i|  >  2|sj. 

Using  this  and  the  estimate  for  F(t],  t)  yields 


|A'(M)I< 

< 


CM 

|s  -  <|  +  1 

CM 

2|s|  +  1  -2  T 


for  all  large  s.  The  estimate  in  the  theorem  follows  easily  by  adjusting  the  constant  as 
necessary. 

For  domains  f!  with  piecewise  C1  boundary,  we  can  show  that  A*  is  The 

proof  is  a  modification  of  VVeyi’s  classical  work  on  the  asymptotics  of  eigenvalues  of 
integral  equations.  The  following  lemma  contains  some  useful  standard  facts  about  the 
Weyl  correspondence.  Again  we  refer  the  reader  to  [8]  for  the  proofs. 

Lemma  6. 

(i)  The  operators  corresponding  to  S(£,  t)  and  S(£  —  —  ^o)  are  unitarily 

equivalent. 

(ii)  Suppose  the  symbols  Si($,t)  and  S2(£,t)  are  related  by  an  orthogo¬ 
nal  change  of  variables.  Then  the  corresponding  operators  are  unitarily 
equivalent. 

We  will  first  prove  that  A*  is  for  domains  ft  of  the  following  form: 


(*)  fi  =  {(£,<)  :<*(/)<£  te[-T,T]} 


6 


where  a(t),/3(t)  are  C1  functions  with  vanishing  endpoint  values.  The  kernel  associated 
to  the  symbol  i(£,0  is: 


’  ir(s  —  t) 

Moreover,  it  is  easy  to  check  that  K  is  a  Lipschitz  function  whose  gradient  DK  exists 
a.e.  and  satisfies  the  inequality 


\DK\  <  C  a.e. 


(2) 


Lemma  7.  If  Q  is  of  the  form  (*)  then  there  are  symmetric  finite  rank  approximations 

ktf(s,t)  =  y^aqXQ 

where  Q  ranges  over  squares  of  the  form 

{(€,0  :  i/N  <£<(i  +  1  )/NJ/N  <t<(j  +  1  )/N)  -  N2  <  i,j  <  N2  -  1 

with  the  property  that 

jj\I<(s,t)-kN(s,t)\2dsdt  is  0(l/N).  (3) 


Proof:  Set 

=w\IIQ 

The  symmetry  of  the  kernel  K  forces  symmetry  of  k^ .  Let  R  =  U<5-  Then 
jjR\K(s,t)-kN(sMdsdt  =  zlL  |A'(s,<)  —  ktf(s,t)\2  dsdt 

-  £  i«i  JJQ  |CA'i2  =  je  JJR 

using  Poincare’s  inequality  [11].  Since  K(s,t)  is  supported  within  the  strip  |s  +  t|  <  2 T, 
equation  (2)  yields  the  estimate 


JJ  \I<(s,t)  -kN(s,t)\2dsdt  < 


C_ 

N' 


(4) 
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The  mean  square  error  over  the  exterior  of  R  can  be  handled  as  follows: 

JljK-k*^JLjKfdsit4  (5) 

in  view  of  lemma  5.  Putting  equations  (4)  and  (5)  together  yields  the  estimate  in 
equation  (3). 

In  view  of  Satz  III  of  Weyl’s  paper  [11],  one  has 

^4W2+i  +  ■•■  <  C/N. 

We  will  apply  this  inequality  to  N/2  where  N  =  \y/k  —  1/2].  (Here  [  ]  denotes  the 
greatest  integer  function.)  This  yields: 

^4(JV/2)2— l  +  '  "  +  ^4^2— l  +  ■  •  •  +  tfc  ■  ■  ■  <  2 C/N. 

Since  the  eigenvalues  have  been  arranged  in  decreasing  order,  the  left-hand  side  of  the 
above  inequality  is  greater  than  or  equal  to  3Ar2A|.  This  yields  an  inequality  of  the 
form  A|  <  C/N3.  Clearly  1/iV  is  0{\/\/k).  These  remarks  imply  that  A*  is  for 

domains  with  property  (*).  We  now  use  a  cut  and  paste  argument  to  derive  the  result 
for  general  Q  with  Cl  boundary. 


Lemma  8.  If  K'{s,t)  and  K"(s,t)  are  two  symmetric  real  kernels  in  L2( R2)  with 

then  the  eigenvalues  of  K(s,t )  =  K'(s,t)  +  K"(s,t)  must  satisfy  the  same  estimate: 
1^*1  ^  pjr t  for  some  constant  C . 


PROOF:  Satz  I  of  Weyl’s  paper  [11]  implies  that 

|A«b|<|Ai|  +  |A2+1| 

IM+il  <  +  lAfc+il 

for  all  positive  integers  k.  From  this  it  is  straightforward  to  verify  that  the  eigenvalues 
of  K(s,t)  have  the  required  decay  property. 

THEOREM  9.  If  Cl  is  a  bounded  domain  with  piecewise  CA  boundary,  then  A*  is 


PROOF:  Clearly,  such  an  fl  can  be  decomposed  into  a  finite  union  of  nonoverlapping 
sub-domains  fl  =  Ujfcfl*  where  each  fl*  can  be  put  into  the  form  (*)  after  a  rigid  motion 
in  the  plane.  By  lemma  6,  each  L\Uk  has  eigenvalues  with  the  sought-after  decay 
property.  Lemma  8  then  implies  that  LXn  also  has  the  same  property. 

This  estimate  is  in  fact  sharp,  at  least  for  annular  regions. 
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Proposition  10.  Let  H  =  {(£,*) :  t  <  (£2  + t2)1!2  <  R}  where  0  <  t  <  R.  Then 

0  <  limsup  &3/4A*  <  oo. 

k— *oo 


PROOF:  According  to  lemma  4,  the  k-th  eigenvalue  is 

rR 


exp(— 27rr2)Lt(47rr2)  rdr. 


The  following  classical  asymptotic  expansion  for  Laguerre  polynomials,  valid  for 
x  €  [e',  i?]  with  0  <  e'  <  R! ,  will  be  essential  [12,  Theorem  8.22.2]: 

jr1/2!1/4^/4  exp(— x/2)  Lk(x)  =  cos(2 (kx)1/2  —  ?r/4)  ^1  +  A\{x)k~1^2  +  Oik-1)  j 

+  sin(2(fca:)1/2  -  tt/4)  (Bi(x)Jr1'2  +  0(fc-1))  • 

Applying  this  formula  with  x  =  47rr2  yields 

|A*|  =  k-^Ioik)  +  k~3'4h(k)  +  0(k~V*)  (6) 


where 


and 


rR 

I0(k)  =  2^^  J  r1!2  cos(4(A:7r)1/2r  —  7r/4)  dr 
rR 

I\(k)  =  23/,27rly/4  /  r1/2 Ai(47rr2)  cos(4(fc7r)1/2r  —  7r/4) dr 

rR 

+  2 3/27r1/4  j  r1//2Bi(47rr2)  sin(4(fc7r)1/,2r  —  7r/4)  dr. 


Consider  the  behavior  of  Io(k).  Using  a  double-angle  formula,  one  has 

R.  R 

/©(*)  =  23/,2it1/4  J  r1^2  cos(4(fc7r)1/,2r)  dr  +  23^27 r1|/4  J  r1^2  sin(4(fc7r)^2r)  dr. 

These  are  the  Fourier  cosine  and  sine  integrals  of  the  smooth  function  r1/2  on  the 
interval  [e,  /?],  evaluated  as  4 {kn)1!2.  As  such,  both  integrals  are  0(k -1/2).  A  similar 
argument  shows  that  the  integrals  in  h{k)  are  0{k -1I2)  as  well. 


9 


Finally,  with  a  little  more  care  one  can  show  that 


limsup &1//2/o(fc)  >  0. 

Jt— oo 

In  fact,  write  the  function  r1/2  =  £(r )+(r^2  — t'(r)),  where  f(r)  is  the  linear  function  that 
interpolates  between  the  endpoint  values  of  r1/2  at  e  and  R.  The  function  r1/2  —  £(r)  is  a 
Lipschitz  function  on  the  real  line  with  support  in  [e,  R).  Consequently,  the  Fourier  sine 
and  cosine  integrals  of  this  function  evaluated  at  4(kir)1/2  are  0(k~l).  It  is  therefore 
enough  to  show  that  the  integral 


A0(k) 


=  3? 


a'*l  4 


rR 

J  ((y)exp(4i(kTr)lt2r) 


has  the  property  that 

limsup  k1^2Ao(k)  >  0. 

it— oc 

This  is  easy  to  show  directly. 

For  general  domains,  we  obtain  the  weaker  result  that  the  sequence  of  eigenvalues 
is  not  absolutely  summable. 

Proposition  11.  The  series  A*  is  not  absolutely  convergent. 

PROOF:  It  is  well  known  that  W(<f>k,<f>i)  is  an  orthonormal  basis  for  L2( R2).  The 
equations  //p  W(<£k,  <£i)  =  ^kl^k  imply  that  xn  =  YlXkW{<t>k,<f>k)  in  L2(R2).  On  the 
other  hand,  if  |A*|  <  oo  then  Yj^kW{(}>k,4>k)  would  have  to  converge  uniformly  to 
an  element  of  Co(R2)  (see  theorem  1  part  2).  This  is  clearly  a  contradiction. 

We  now  examine  the  smoothness  and  decay  of  the  eigenfunctions.  We  assume  that 
f]  is  an  open  set  contained  in  a  rectangle  \—B,B]  x  [— T,  T]  of  which  all  cross-sections 
in  the  £  and  t  directions  consist  of  at  most  M  intervals.  We  now  examine  an  equation 
of  the  form  LXn<f>  =  \<j>,  where  A  is  a  nonzero  real  number  and  <f>  G  T2(R).  Fubini’s 
theorem  implies  that 

A ^(s)  =  J  K(s,t)<f>(t)  dt 

holds  for  all  s  €  R\Z,  where  Z  is  some  set  of  measure  zero.  Lemma  5  yields  the  estimate 


< 

< 


c 

i  +  M 

CyfiT 

M  +  1 


/: 


-s+2T 
2  T 


W&Wl2  < 


\m\dt 

C' 

W  +  i 
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for  all  s  6  R\Z.  This  last  estimate  implies  that  <f>  €  L°°( R).  Now  define 

£(s)  =  sup  \<f>{s)\. 

M>W 

Note  that  £(s)  is  even  and  decreasing  in  )sj.  The  quantity  )A<^(s)|  can  now  be  estimated 
as  follows: 

r <  f—s+2T 

|A^)|  < 


< 


w  + 

ACT 


H  +  1 


/-s+il 

\4>{t)\dt 

s-2T 

£(\s\-2T) 


for  all  s  e  (R\-Z)  D  {.s'  :  js'j  >  2 T}.  Therefore,  given  any  b  >  1,  there  is  an  .so  such  that 
if  \s\  >  |s0|  and  s  €  R\Z  then 


Iterating  this  estimate  yields 


€(\s\  +  2T)  <  ^(|s|). 


e(\s\  +  2nT)<-£(\s\). 


(7) 


It  is  straightforward  to  check  that 

|^(s)|  <  C6“|s|/(2T)  Vs  6  R\Z. 

As  a  consequence,  <)>  is  smooth  and  has  all  its  derivatives  in  L2 . 

Now,  by  theorem  1  part  4,  we  have  that 

ffmh,4>i)=  IS  W(tk,4>i)  =  Wki 

J Ju  J Jn 

where  U  =  {(<^,  <)  :  (— 1,£)  €  I)}.  Hence  Lx^<j>k  =  A*<^t,  and  by  the  preceding  discussion 
4>  is  smooth  and  has  all  its  derivatives  in  Ir .  Now  for  a  given  sj, 

=  f 

3=z$\  * 

Using  the  estimate  analogous  to  equation  (7)  for  4>  with  6  >>  1,  one  has  |<^(£)|  < 
Ce~u\t I.  Then 

Jo 


n!  dsn 


s=sl 


< 


cr(n) 


-  n! 2n+1  2"+1n 

Hence  the  power  series  of  <f>  at  s  =  sj  converges  in  some  interval  around  sj.  Because  of 
the  symmetry  in  the  role  of  <f>  and  0,  the  same  observations  hold  for  <j>.  The  preceding 
discussion  is  summarized  in  the  following  theorem. 
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THEOREM  12.  Suppose  fl  is  an  open  set  contained  in  the  interval  [— B,  B]  x  [— T,  T\ 
with  the  property  that  all  its  cross-sections  in  both  the  £  and  t  directions  consist  of  at 
most  M  intervals.  Then: 

(i)  for  any  b  >  0,  there  is  a  constant  Cj  such  that 

|<M-S)|  5:  Vs,  and 

imi<C>e-M  V£, 

and 

(ii)  <f>  and  <p  are  analytic  and  have  all  their  derivatives  in  L2 . 

Note  that  (1)  actually  implies  (2)  in  theorem  12  by  the  Paley- Wiener  theorems 
[13],  although  we  have  chosen  to  give  the  elementary  argument. 

It  is  instructive  to  consider  the  case  when  H  is  a  ball  centered  at  the  origin.  It  is 
well  known  that,  in  this  case,  the  are  Hermite  functions  (see  lemma  4).  The  Hermite 
functions  hj(t)  in  lemma  4  will  satisfy  estimates  of  the  form  hj(t)  <  Ce~^~e^  for 
any  e  >  0.  It  is  directly  evident  that  they  will  then  satisfy  the  weaker  inequalities  in 
theorem  12.  This  theorem  states  that  this  weaker  decay  statement  holds  for  general 
regions  in  the  time-frequency  plane.  We  do  not  know  whether  these  estimates  can  be 
improved. 
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SECTION  4 

NUMERICAL  EXAMPLES 


For  illustration  purposes,  we  provide  four  numerical  examples  of  time-frequency 
localization.  These  examples  are  obtained  by  discretizing  the  kernel  in  the  integral 
representation  of  the  operator  given  in  lemma  5.  We  consider  localization  on  domains 
of  the  form  of  a  “zigzag,”  a  disk,  a  rectangle,  and  a  parallelogram,  as  indicated  in 
figure  1.  Note  that  while  the  operators  here  are  all  Hermitian,  so  that  they  have  real 
eigenvalues,  they  do  not  generally  have  real  eigenfunctions  unless  their  kernels  are  real. 
This  is  true  if  the  symbol  S(£,  t)  satisfies  S(—£,t)  =  S((,t).  For  clarity  of  eigenfunction 
representation,  projection  domains  were  chosen  to  allow  this  symmetry  when  possible. 

Localization  on  the  unusual  zigzag  domain  (figure  la)  produces  the  unfamiliar 
eigenfunctions  in  figure  2  (a s  these  are  complex,  we  have  plotted  their  magnitude).  Eight 
samples  per  unit  length  are  used,  over  the  time  domain  [—8,8].  Most  of  the  energy  is 
concentrated  in  the  desired  region  [—2,2],  although  there  is  some  leakage.  Note  that  the 
first  three  eigenfunctions  have  one,  two,  and  three  peaks,  respectively.  For  the  case  of 
the  disk  (figure  lb),  a  plot  of  the  eigenvalues  is  shown  in  figure  3.  As  noted  in  lemma  4, 
the  eigenfunctions  in  this  case  are  the  well-known  hermite  functions.  Next,  localization 
on  the  rectangle  (figure  lc)  allows  comparison  to  the  prolate  spheroidal  wavefunctions. 
The  first  nine  solutions  for  the  prolate  spheroidal  and  Weyl  operator  cases  are  depicted 
in  figures  4  and  5.  Again,  eight  samples  per  unit  length  are  used  over  the  time  domain 
[—8,8].  While  there  is  energy  leakage  outside  the  desired  domain  [—2,2],  the  amount 
differs  in  the  two  cases.  The  prolate  spheroidal  wavefunctions  have  \/x  decay  while 
the  Weyl  eigenfunctions  have  exponential  decay.  This  difference  is  visible  when  one 
compares  the  last  three  eigenfunctions  in  each  case.  A  plot  of  the  eigenvalues  is  also 
included. 

Finally,  we  provide  a  simple  illustration  of  these  ideas  in  the  context  of  filtering 
Gaussian  noise  from  a  corrupted  linear  FM  (chirp)  signal.  Figures  6a  and  7a  show, 
respectively,  the  Wigner  distribution  intensity  plot  and  the  actual  plot  of  the  real  part 
of  a  linear  FM  chirp  centered  at  0  frequency.  (Although  chirp  signals  used  in  radar  are 
not  centered  at  0  frequency,  that  is  irrelevant  for  our  purposes  because  of  the  covariance 
of  the  Wigner  distribution  under  time  and  frequency  shifts  (theorem  1,  part  5).)  In 
the  discretization,  16  samples  per  unit  length  are  used  over  the  time  domain  [—2,2]. 
No  windowing  is  applied  to  remove  the  echo  effect  in  these  Wigner  plots,  although  we 
have  found  a  simple  cosine-squared  window  to  be  effective.  Figures  6b  and  7b  show 
this  signal  after  0  dB  Gaussian  white  noise  has  been  added.  To  filter  the  noise,  we  note 
that  theoretically,  the  Wigner  distribution  of  a  chirp  signal  is  a  measure  concentrated 
along  a  diagonal  line  corresponding  to  the  slope  of  the  chirp.  In  particular,  a  chirp  can 
be  localized  in  any  domain  containing  its  time- frequency  support. 
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Figure  1.  Localization  Domains:  a)  Zigzag  Domain, 
b)  Disk  of  Area  50,  c)  4x2  Rectangle,  and 
d)  Parallelogram 


As  an  elementary  example  of  time-frequency  localization,  the  noisy  chirp  of  figure 
7b  is  projected  onto  the  first  two  eigenfunctions  (weighted  according  to  the  eigenvalues) 
for  the  domain  in  figure  Id.  This  results  in  the  signal  in  figure  7c,  with  a  Wigner 
distribution  as  shown  in  figure  6c.  A  plot  of  the  eigenvalues  for  this  localization  operator 
is  provided  in  figure  7d.  In  fact,  the  chirp  in  figure  7a  is  orthogonal  to  the  second 
eigenfunction,  illustrating  an  interesting  fact.  Numerically,  the  chirp  appears  to  be  an 
exact  solution  to  the  problem  of  localizing  onto  an  infinite  diagonal  band  domain.  For 
more  numerical  examples  of  time-frequency  localization,  the  reader  is  referred  to  [6,7]. 
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Figure  2.  Magnitude  of  Eigenfunctions  1-3  for  Localization 
onto  a  Zigzag  Domain  (in  figure  la),  with  a  Plot  of  Eigenvalues 
Sorted  According  to  Decreasing  Magnitude 

It  is  interesting  at  this  point  to  remark  on  the  conjecture  of  Flandrin  [5]  that 
for  localization  onto  convex  domains,  the  top  eigenvalue  is  bounded  above  by  1.  This 
seems  to  hold  (at  least  numerically)  in  our  examples,  even  for  figure  la,  which  is  not 
convex.  However,  something  like  convexity  is  certainly  necessary  in  general,  since  we 
have  nonconvex  numerical  examples  where  the  top  eigenvalue  exceeds  4. 
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Figure  3.  Weyl  Eigenvalue  Plot  for  Localization 
onto  a  Disk  of  Area  50  (figure  lb),  Sorted 
According  to  Decreasing  Magnitude 
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4.  Prolate  Spheroidal  Wavefunctions  1-9  for  Localization 
on  a  Rectangle  (figure  lc),  with  Eigenvalue  Plot 
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Figure  5.  Weyl  Eigenfunctions  1-9  for  Localization 
on  a  Rectangle  (figure  lc),  and  a  Plot  of 
Eigenvalues  in  Decreasing  Magnitude 
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Frequency 


Frequency 


Figure  6.  Wigner  Distribution  Intensity  Plot 
of  (a)  Chirp  Signal,  (b)  Chirp  in  0  dB  Gaussian 
White  Noise,  and  (c)  Filtered  Noisy  Signal 


(c)  (d) 

Figure  7.  Plots  of  the  Real  Part  of  (a)  Chirp  Signal, 
(b)  Chirp  Signal  in  0  dB  Gaussian  White  Noise,  and 
(c)  Filtered  Noisy  Signal;  Part  (d)  Is  an  Eigenvalue  Plot. 
The  First  Two  Eigenfunctions,  Weighted  by  their 
Eigenvalues,  Are  Used  in  the  Filtering. 
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